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Synthesis  of  the  spinel  lithium  titanate  Li4Ti50i2  by  an  alkoxide-free  sol-gel  method  is  described.  This 
method  yields  highly  pure  and  crystalline  Li4Ti50i2  samples  at  relatively  low  temperature  (850 °C)  and 
via  short  thermal  treatment  (2h).  6Li  magic  angle  spinning  nuclear  magnetic  resonance  (MAS  NMR) 
measurements  on  these  samples  were  carried  out  at  high  magnetic  field  (21.1  T)  and  over  a  wide  tem¬ 
perature  range  (295-680  K).  The  temperature  dependence  of  the  chemical  shifts  and  integral  intensities 
of  the  three  6Li  resonances  demonstrates  the  migration  of  lithium  ions  from  the  tetrahedral  8 a  to  the 
octahedral  16c  sites  and  the  progressive  phase  transition  from  a  spinel  to  a  defective  NaCl-type  struc¬ 
ture.  This  defective  structure  has  an  increased  number  of  vacancies  at  the  8a  site,  which  facilitate  lithium 
diffusion  through  16c^  8a  ^  16c  pathways,  hence  providing  an  explanation  for  the  reported  increase 
in  conductivity  at  high  temperatures.  Molecular  dynamics  simulations  of  the  spinel  oxides  Li4+xTi50i2, 
with  0<x<3,  were  also  performed  with  a  potential  shell  model  in  the  temperature  range  300-700 K. 
The  simulations  support  the  conclusions  drawn  from  the  NMR  measurements  and  show  a  significant 
timescale  separation  between  lithium  diffusion  through  8a  and  16c  sites  and  that  out  of  the  16d  sites. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Titanium  oxide  polymorphs  have  received  much  attention  due 
in  part  to  their  application  in  solid-state  lithium-ion  batteries 
[  1  -3 ].  In  particular,  spinel-structured  lithium  titanates  have  shown 
promising  properties  for  use  as  electrode  materials  [4,5].  For  exam¬ 
ple,  the  spinel  Li4Ti50i2  (Fd3m)  is  being  considered  as  a  potential 
anode  material  for  commercial  purposes  due  to  its  excellent  elec¬ 
trochemical  properties,  which  include  a  low  reduction  potential 
(~1.5  V  vs.  Li/Li+),  high  reversibility  of  the  Li+  insertion/extraction 
process,  and  a  flat  operating  voltage  [6].  In  addition,  this  mate¬ 
rial  undergoes  minimal  structural  changes  upon  lithium  insertion 
and  extraction  and  has  been  termed  a  zero-strain  insertion  mate¬ 
rial  [6,7].  Indeed,  in-situ  energy  dispersive  X-ray  diffraction  shows 
that  variations  in  lattice  parameters  upon  lithium  insertion  and 
extraction  are  at  most  1%  [8].  Finally,  Li4Ti50i2  can  be  charged 
up  to  Li7Ti50i2,  which  corresponds  to  a  theoretical  capacity  of 
175mAhg_1,  and  shows  an  excellent  capacity  maintenance  ratio 
of  99%  even  after  the  100th  cycle  [9].  It  should  be  noted  that  it  has 
been  suggested  based  on  both  experimental  [10,11]  and  computa¬ 
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tional  [12]  results  that  further  lithium  can  be  inserted  in  Li7Ti50i2 
up  to  a  stoichiometry  of  LigTi50i2  at  a  potential  of  less  than  1 V, 
resulting  in  a  theoretical  capacity  of  293  mA  h  g-1 .  These  excellent 
electrochemical  properties  make  this  material  a  prime  candidate 
for  use  as  anode  material  in  commercial  lithium-ion  batteries  [13]. 

The  use  of  Li4Ti50i2  in  commercial  devices  will  be  facilitated 
by  a  synthesis  method  that  can  easily  be  implemented  for  bulk 
production.  Flowever,  the  synthesis  techniques  reported  to  date, 
such  as  solid-state  synthesis  [6,14,15],  sol-gel  [16,17],  combus¬ 
tion  technique  [18,19],  rheological  phase  reaction  [20]  and  ball 
milling  method  [15,21],  for  preparing  Li4Ti50i2  use  excessive  heat 
treatment  (>850  °C  and  >1 2  h)  and/or  highly  toxic  and  air  sensitive 
starting  materials,  which  complicates  their  use  in  industrial  envi¬ 
ronments.  In  this  report,  we  present  a  simple  sol-gel  method  using 
water  and  a  soluble  titanium  citrate  complex  as  starting  material, 
which  helps  to  prepare  Li4Ti50i2  in  open  air  and  relatively  short 
heat  treatment  (2  h). 

The  attractive  electrochemical  properties  of  Li4Ti50i2  arise 
mainly  from  the  lithium  diffusion  behavior  in  its  crystal  lattice.  Con¬ 
siderable  efforts  have  been  made  to  understand  lithium  diffusion 
and  local  structure  in  Li4Ti50i2.  The  crystal  structure  of  Li4Ti50i2 
was  initially  reported  by  Deschanvres  et  al.  [22]  and  subsequently 
by  Kataoka  et  al.  [23].  The  latter  structure  is  shown  in  Fig.  1.  Note 
that  in  Fig.  1  all  possible  lithium  sites  are  shown  as  occupied  for  pre¬ 
sentation  purposes.  The  distribution  of  lithium  ions  between  these 
sites  is  discussed  in  what  follows.  In  Li4Ti50i2,  the  octahedral  16d 
sites  are  randomly  occupied  by  lithium  and  titanium,  the  tetra- 
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Fig.  1.  Top  figure:  crystal  structure  of  LUTisO^  showing  the  three  lithium  sites: 
tetrahedral  8a  (white),  octahedral  16c  (grey)  and  octahedral  16d  (black).  Tita¬ 
nium  and  oxygen  atoms  are  shown  as  small  dark  grey  and  large  light  grey  balls, 
respectively.  All  three  sites  are  shown  as  occupied  for  presentation  purposes.  The 
occupancy  of  these  sites  will  vary  with  Li  content  and  temperature,  as  discussed  in 
the  text.  Bottom  figure:  inter-atomic  distances  between  16d-16c  (2.95  A),  16c-8a 
(1.81  A),  and  8a-16d  (3.46  A)  site  pairs  based  on  the  structure  reported  by  Kataoka 
et  al.  [23]. 


charge  potential,  whereby  Li4Ti50i2  and  Li7Ti50i2  domains  coex¬ 
ist  and  the  reaction  front  advances  in  the  material  as  the  reaction 
proceeds  [5].  Scharner  et  al.  [24]  provided  evidence  of  a  two-phase 
transformation  upon  lithium  insertion  based  on  X-ray  diffraction 
(XRD)  data  and  analysis  of  high-angle  reflections,  which  allowed 
for  the  determination  of  a  linear  variation  of  the  proportion  of  the 
two  phases  with  increasing  lithium  content.  Interestingly,  Wage- 
maker  et  al.  have  shown  using  neutron  and  X-ray  diffraction  that  in 
the  absence  of  an  applied  potential  the  two-phase  structure  relaxes 
to  a  solid  solution  [25].  This  phenomenon  was  previously  alluded 
to  by  Ronci  et  al.  [26]  on  the  basis  of  changes  in  lattice  parameters 
determined  from  in-situ  synchrotron  radiation  XRD.  A  later  study 
by  Wagemaker  et  al.  revealed  8a  (Li4Ti50i2)  and  16c  (Li7Ti50i2) 
domain  sizes  on  the  order  of  a  few  nanometers  [27]. 

Li4Ti50i2  also  shows  interesting  properties  at  high  tem¬ 
peratures  [28,29].  For  example,  Li4Ti50i2  undergoes  a  nearly 
six-order-of-magnitude  increase  in  conductivity  (~10-3Scm-1) 
at  high  temperatures  (>600  K)  compared  to  the  low  conductivity 
(~10-8  S  cm-1 )  observed  at  room  temperature  [29].  This  enormous 
increase  in  conductivity  has  been  interpreted  as  the  formation 
of  a  disordered  state  at  high  temperatures,  whereby  lithium  ran¬ 
domly  distributes  between  the  8 a  and  16c  sites.  The  displacement 
of  lithium  ions  from  the  8 a  to  the  16c  sites  with  temperature  and 
the  formation  of  a  disordered  state  have  been  inferred  on  the  basis 
of  changes  in  the  intensities  of  infrared  and  Raman  modes  with 
temperature  [28,29]  but  has  not  yet  been  confirmed  with  other 
experimental  probes. 

To  provide  further  insights  into  the  changes  in  lithium  site 
distribution  and  diffusion  behavior  in  Li4Ti50i2  with  tempera¬ 
ture,  we  have  carried  out  6Li  magic  angle  spinning  (MAS)  nuclear 
magnetic  resonance  (NMR)  measurements  in  combination  with 
molecular  dynamics  (MD)  simulations.  Li4Ti50i2  has  previously 
been  studied  with  7Li  NMR;  however,  the  data  in  these  studies 
were  acquired  at  low  magnetic  fields  (<11.7T),  and  hence  might 
have  suffered  from  poor  spectral  resolution  [30,31].  Furthermore, 
these  measurements  were  carried  out  at  temperatures  of  at  most 
500  K  and  were  therefore  not  able  to  probe  the  disordered  phase 
and  related  increase  in  conductivity  above  600  K.  In  the  present 
work,  variable-temperature  6Li  MAS  NMR  measurements  were 
performed  up  to  700  K  and  under  high  magnetic  field  (21.1  T)  to 
explore  the  lithium  local  structure  and  diffusion  in  Li4Ti50i2  at 
high  temperatures.  Although  NMR  can  offer  a  conspicuous  view 
of  the  diffusion  properties  of  complex  materials,  interpretation  of 
NMR  data  is  theoretically  complex  and  thus  often  relies  on  com¬ 
plementary  information  from  computational  models  to  elucidate 
the  ion  distribution  and  dynamics.  For  example,  in  a  recent  study 
[32],  results  obtained  from  MD  simulations  helped  assign  the  three 
resonance  peaks  obtained  in  NMR  measurements  of  L^TiOs.  In 
this  work,  experimental  NMR  data  is  combined  with  results  from 
molecular  dynamics  simulations  to  determine  the  lithium  distri¬ 
bution  and  diffusion  in  Li4Ti50i2  at  high  temperatures. 


hedral  8a  sites  are  occupied  by  lithium  only  and  the  16c  sites  are 
vacant.  Upon  lithium  insertion,  the  lithium  ions  move  from  the  8a 
to  the  16c  sites  and,  at  a  composition  of  Li7Ti50i2,  a  rocksalt  struc¬ 
ture  is  formed  whereby  all  16c  sites  are  occupied  and  all  8a  sites 
are  empty.  Lithium  insertion  is  therefore  proposed  to  be  limited  by 
the  number  of  available  16c  sites.  Flowever,  as  already  mentioned, 
higher  lithium  contents  have  been  reported  [10-12],  whereby  the 
8 a  sites  begin  to  fill  once  all  1 6c  sites  are  occupied.  Such  an  insertion 
behavior  suggests  that  lithium  insertion  is  limited  by  the  number 
of  tetravalent  titanium  ions  rather  than  by  the  number  of  available 
octahedral  sites. 

Li4Ti50i2  shows  a  remarkably  flat  charge  and  discharge  poten¬ 
tial  of  approximately  1 .55  V  [4-6].  Therefore,  lithium  insertion  was 
suggested  to  follow  a  two-phase  reaction  based  on  the  observed  flat 


2.  Experimental 

The  spinel  type  Li4Ti50i2  material  is  prepared  by  modified 
pechini  type  sol-gel  method  using  non-alkoxide  materials  simi¬ 
lar  to  previous  report  [33].  To  prepare  1  g  of  Li4Ti50i2,  0.39 g  of 
Ti  metal  powder  was  dissolved  in  an  aqueous  solution  contain¬ 
ing  20  mL  of  H202  (30%)  and  5mL  of  ammonia  (35%)  in  an  ice 
bath.  This  yielded  a  yellowish  transparent  titanium  peroxo  solution. 
Citric  acid  (15  g)  was  then  added  in  excess  to  the  titanium  per¬ 
oxo  solution  to  form  a  metal  chelate.  Fligh-purity  Li2C03  (0.246  g) 
was  dissolved  in  denaturated  ethyl  alcohol.  This  solution  was  then 
added  to  the  titanium  peroxo  solution  and  stirred  well  to  ensure  a 
homogeneous  distribution.  Ethylene  glycol  (20  g)  was  then  added 
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Fig.  2.  XRD  patterns  of  powder  precursor  samples  heat-treated  at  different  tem¬ 
peratures  (a)  600 °C  for  2  h,  (b)  700 °C  for  2  h,  (c)  850  °C  for  2  h.  The  spinel  Li4Ti50i2 
diffraction  lines  are  matched  with JCPDS  data  (#  049-0207)  and  the  rutile  Ti02  impu¬ 
rity  phase  indicated  by  the  asterisk  (*)  are  identified  using  JCPDS  data  (#021-1276). 


to  the  above  solution,  which  became  reddish  in  color.  The  tempera¬ 
ture  was  then  increased  to  150  °C  for  polymerization,  which  yields 
chemically  stable  black  polymer  precursor.  Therefore,  pyrolysis  of 
the  polymer  was  performed  at  350  °C  that  yielded  a  black  powder 
precursor.  This  final  powder  precursor  was  subjected  to  calcina¬ 
tions  in  open  air  for  2  h  at  temperatures  between  600  and  850  °C. 
The  resulting  product  from  this  sol-gel  method  was  analyzed  by 
X-ray  powder  diffractometry  using  a  Philips  Xpert  X-ray  diffrac¬ 
tometer  using  CuKa  radiation  at  X  ~  1.54  A.  Fig.  2  shows  the  XRD 
pattern  of  the  precursor  after  heat-treatment  at  different  temper¬ 
atures.  The  powder  precursor,  when  heat-treated  below  600  °C,  is 
amorphous  in  nature.  Its  gray  color  indicates  the  presence  of  resid¬ 
ual  carbon.  When  heat-treated  at  700  °C  the  powder  appears  white 
in  color,  which  suggests  the  complete  burnout  of  residual  carbon 
from  the  precursor  powders.  The  formation  of  the  crystalline  spinel 
phase  (JCPDS#  049-0207)  along  with  rutile  (JCPDS#021-1276)  was 
observed  in  the  XRD  patterns  (see  Fig.  2)  for  samples  heat-treated 
at  700  °C.  The  highly  crystalline  and  pure  spinel  phase  has  been 
identified  for  the  samples  heat-treated  at  850  °C  for  2  h.  The  crucial 
step  adopted  in  our  present  study  was  to  prepare  a  highly  water- 
soluble  starting  precursor  to  avoid  the  use  of  alkoxides  reported 
in  previous  literatures  [17,18,34]  which  are  not  stable  in  open  air. 
This  was  achieved  by  dissolving  Ti  metal  in  a  solution  of  hydrogen 
peroxide  and  ammonia,  as  previously  reported  [35]  which  yield  an 
air-stable  and  water-soluble  peroxo-citrato-titanium  complex  as 
titanium  source  for  sol-gel  chemistry.  Additionally,  the  relatively 
low  temperature  and  the  short  time  used  in  this  synthesis  minimize 
the  lithium  loss,  as  it  is  generally  observed  in  the  high-temperature 
heat  treatments  in  other  synthesis  procedures  [6,1 9].  Therefore,  the 
modified  sol-gel  method  used  in  this  study  offers  an  inexpensive 
and  convenient  aqueous  synthesis  approach  that  is  scalable  to  meet 
the  need  of  large  scale  commercial  applications. 

6Li  NMR  measurements  were  performed  on  a  Varian  900  Inova 
magnet  (B0  =  21.1  T  and  6Li  Larmor  frequency  of  132.56  MFIz)  with 
magic  angle  spinning  (MAS)  at  7kFIz  using  5  mm  zirconia  rotors 
and  a  home-made  high-temperature  MAS  probe.  The  6Li  spectra 
are  recorded  with  single  pulse  measurements  and  the  necessary 
recycle  delay  (di )  was  identified  using  array  of  d\  varying  from  1  to 
700  s  at  each  temperature.  For  the  temperature  range  of  our  study, 


500  s  recycle  delay  was  sufficient  to  acquire  quantitative  spectra 
at  all  temperatures.  Chemical  shifts  were  externally  referenced  to 
a  1  M  LiCl  solution  (5iso  =  0  ppm).  The  estimated  uncertainties  in 
chemical  shifts,  determined  through  calibration  using  the  same 
reference  sample,  were  0.1  ppm.  Accurate  chemical  shifts  and  inte¬ 
gral  intensities  of  different  contributions  in  each  spectrum  were 
obtained  from  fitting  the  line  shapes  of  the  resonance  lines  with  a 
mix  of  Gaussian  and  Lorentzian  line-shapes  using  the  DMFIT  pro¬ 
gram  [36].  All  the  6Li  NMR  measurements  were  carried  out  on 
powder  Li4Ti50i2  sample  without  any  6Li  isotope  enrichment. 

3.  Computational  methods 

3.1.  Potential  model 


The  calculations  reported  in  this  work  are  based  on  the  Born 
model  of  solids  [37],  in  which  the  atoms  of  a  system  are  represented 
as  point-charge  particles  that  interact  via  long-range  Coulom- 
bic  forces  and  short-range  interactions.  The  latter  are  described 
by  parameterized  functions  and  represent  the  repulsion  between 
electron-charge  clouds  and  the  van  der  Waals  attraction  forces.  In 
this  work,  the  short-range  interactions  are  described  using  a  Buck¬ 
ingham  potential  and  therefore  the  pair-wise  interaction  energy 
takes  the  following  form: 


UU  = 


1  M] 
4tT£0  Ty 


+  Ay  exp 


(1) 


where  the  first  term  is  the  electrostatic  potential  energy  and 
the  second  and  third  terms  represent  the  electron-charge  cloud 
repulsion  and  van  der  Waals  attraction  energies,  respectively.  The 
parameters  (q,  A,  p,  and  C)  used  in  this  study  for  titanium  and  oxy¬ 
gen  ions  are  those  optimized  by  Matsui  and  Akaogi  (MA)  [38  ]  for  the 
Ti02  polymorphs  rutile,  anatase,  brookite,  and  Ti02-II.  The  poten¬ 
tial  parameters  used  to  describe  lithium  ions  were  determined 
by  fitting  the  potential  parameters  for  Li  to  the  lattice  parame¬ 
ters,  lattice  constants,  and  bulk  modulus  of  Li20  while  keeping  the 
MA  potential  parameters  constant.  The  lithium  partial  charge  was 
determined  from  the  Li20  stoichiometry  while  keeping  the  oxygen 
partial  charge  from  the  MA  model  constant  (i.e.,  qu  =  -q0/2  =  +0.549 
e).  Details  of  the  derivation  of  the  Li  parameters  are  provided  in 
Supplemental  information.  In  addition,  our  model  accounts  for  the 
polarizability  of  oxygen  ions  via  a  shell  model  first  introduced  by 
Dick  and  Overhauser  [39].  In  the  shell  model,  a  polarizable  ion  is 
composed  of  two  particles,  a  core  and  a  shell,  which  share  the  ion’s 
charge  and  are  linked  by  a  harmonic  spring,  k: 


Uc-s  =  kx  r(Ls 


(2) 


where  rc_s  is  the  core-shell  separation  distance.  The  derivation  of 
the  shell  model  version  of  the  MA  potential  is  described  in  a  pre¬ 
vious  publication  [40].  This  shell  model  has  proved  to  be  a  reliable 
and  transferable  potential  model.  For  example,  it  gives  an  excellent 
description  of  the  crystal  structure  of  several  Ti02  polymorphs  and 
reproduces  lattice  distortions,  obtained  from  density  functional 
theory  calculations,  due  to  the  introduction  of  hole  and  electron 
polarons  in  rutile  and  anatase  crystal  lattices  [40].  In  addition,  the 
potential  model  yielded  excellent  agreement  with  the  experimen¬ 
tal  crystal  structure  of  Li2Ti03  and  helped  provide  an  explanation 
for  the  presence  of  a  NMR  peak  with  a  chemical  shift  correspond¬ 
ing  to  a  lithium  tetrahedral  site  whereas  the  pure  crystal  structure 
only  contains  lithium  ions  in  octahedral  coordination  [32].  All  the 
potential  parameters  used  in  this  work  are  summarized  in  Table  1. 

The  ability  of  the  potential  model  for  simulating  Li4Ti50i2  was 
evaluated  by  comparing  the  calculated  and  experimental  lattice 
parameters.  A  200-ps  NPT  molecular  dynamics  simulation  at  300  K 
with  the  computer  program  DL_POLY  [41]  yielded  a  =  8.443  A.  X- 
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Table  1 

Potential  parameters  used  in  this  work.  The  oxygen  core-shell  spring  constant,  k,  is 
set  to  44.3  eVA-2. 


Core  and  Shell  chargesa 

Species 

Core (e) 

Shell  (e) 

Ti4+ 

2.196 

_ 

Ti3+ 

1.647 

- 

Li+ 

0.549 

- 

o2- 

0.500 

-1.598 

Buckingham  potential  parameters b 

Ion  pair  (y) 

Ay  (eV) 

Pij  (A) 

CtJ  (evA6) 

Li-Li 

38533.955 

0.100 

0.00 

Ti-Li 

33089.570 

0.127 

0.00 

Ti-Ti 

31120.528 

0.154 

5.25 

Ti-0 

16957.710 

0.194 

12.59 

Li-0 

15465.549 

0.167 

0.00 

0-0 

11782.884 

0.234 

30.22 

a  Core-shell  interaction  potential  form:  V  =  k  x  r2_s. 
b  Buckingham  potential  form:  Vy  =  Ay  x  exp(—ry/py)  -  Cy-r6. 


ray  diffraction  data  of  Deschanvres  et  al.  [22]  and  Kataoka  et  al.  [42] 
gave  a  =  8.357  A  and  a  =  8.352  A,  respectively.  Wagemaker  et  al.  [25] 
reported  a  value  of  the  lattice  parameter  of  8.3595  A  from  neutron 
diffraction  data.  The  calculated  lattice  parameter  is  in  good  agree¬ 
ment  with  experimental  data  with  an  error  of  at  most  1.1%.  The 
potential  model  predicts  a  linear  decrease  in  a  with  increasing  Li 
content  for  0  <x<3  in  Li4+xTi50i2  in  agreement  with  experimen¬ 
tal  data  [25].  For  Li7Ti50i2,  the  molecular  dynamics  simulations 
predict  a  =  8.41 1  A  or  a  reduction  of  0.4%,  which  is  larger  than  that 
observed  by  Wagemaker  et  al.  [25  ]  (just  under  0.1  %)  but  of  the  same 
order  of  magnitude  and  similar  (0.7%)  to  that  calculated  by  Zhong 
et  al.  [12]  using  planewave  density  functional  theory. 

3.2.  Molecular  dynamics  simulations 

MD  simulations  were  carried  out  with  the  computer  program 
DL_POLY  [41  ]  at  zero-applied  pressure.  Equilibration  runs  were  ini¬ 
tially  performed  in  the  NPT  ensemble  (constant  number  of  particles, 
constant  pressure,  and  constant  temperature),  whereby  both  the 
cell  shape  and  size  were  able  to  vary,  for  1  Ops  to  obtain  the  aver¬ 
age  cell  dimensions.  Production  runs  were  then  carried  out  in  the 
NVT  ensemble  (constant  number  of  particles,  constant  volume,  and 
constant  temperature)  for  10  ns  with  the  cell  dimensions  fixed  to 
the  average  values  obtained  in  the  equilibration  runs.  This  was  done 
to  facilitate  the  mean  square  displacement  analysis.  In  these  sim¬ 
ulations,  the  temperature  and  pressure  were  kept  constant  by  use 
of  the  Nose-Hoover  thermostat  [43]  and  the  Hoover  barostat  [44], 
respectively.  The  electrostatic  forces  were  calculated  by  means  of 
the  Ewald  summation  method  [45].  A  9  A  cutoff  was  used  for  the 
short-range  interactions  and  the  real  part  of  the  Ewald  sum.  The 
Verlet  leapfrog  integration  algorithm  was  used  to  integrate  the 
equations  of  motion  with  a  time  step  of  0.2  fs.  The  shells  were  given 
a  mass  of  0.2  a.u.  and  their  motion  treated  as  that  of  the  cores  fol¬ 
lowing  the  adiabatic  shell  model  first  introduced  by  Mitchell  and 
Fincham  [46]. 

4.  Results  and  discussion 

4.1.  High-temperature  6 Li  MAS  NMR 

Fig.  3  shows  the  temperature  dependence  of  6Li  MAS  NMR  spec¬ 
tra  ofLi4Ti5012  acquired  at  a  magnetic  field  of  21 .1 T  using  a  sample 
spinning  speed  of  7  kHz.  At  295  K  (bottom  spectrum  in  Fig.  3),  the 
Li4Ti50i2  sample  shows  three  overlapping  resonances  at  0.65, 0.1 7, 
and  -0.45  ppm  indicating  three  different  crystallographic  posi- 


Fig.  3.  (a)  6Li  MAS  NMR  spectra  measured  with  a  21.1  T  magnetic  field  and  spin¬ 
ning  speed  of  7  kHz  for  Li4Ti50i2.  The  peaks  are  labeled  with  Wyckoffs  notation  for 
different  lithium  positions  in  Li4Ti50i2. 


tions  for  lithium  in  Li4Ti50i2,  which  agrees  with  single-crystal  data 
[23,42].  Fig.  1  shows  the  crystal  structure  of  Li4Ti50i2  as  reported 
by  Kataoka  et  al.  [23,42]  with  all  possible  lithium  sites  occupied. 
As  previously  explained,  Li4Ti50i2  has  two  octahedral  and  one 
tetrahedral  lithium  positions,  which  are  labeled,  in  Wyckoffs  nota¬ 
tion,  16c,  16d,  and  8a,  respectively.  The  6Li  NMR  peak  assignments 
were  made  based  on  previous  reports  of  lithium  insertion  in  tita¬ 
nium  oxide  materials  [47-50]  and  6Li  chemical  shift  correlation 
data  [51  ].  According  to  6Li  chemical  shift  correlation  data  [51  ],  the 
nature  of  the  lithium  coordination  environment  is  the  main  driving 
force  for  the  observed  chemical  shift.  For  example,  the  chemi¬ 
cal  shift  (with  respect  to  1  M  aqueous  LiCl  solution)  of  lithium  in 
tetrahedral  coordination  is  expected  around  1  ppm  towards  higher 
frequencies  (positive  shift)  whereas  that  of  lithium  in  octahedral 
coordination  is  expected  around  -1  ppm  (negative  shift).  Follow¬ 
ing  this  chemical  shift  correlation  data,  the  tetrahedral  position  (8a) 
can  be  assigned  to  the  6Li  resonance  at  0.65  ppm  whereas  the  two 
octahedral  lithium  positions,  i.e.,  16d  and  16c,  can  be  assigned  to 
the  6Li  resonances  observed  at  0.17  and  -0.45  ppm,  respectively, 
based  on  previous  reports  [47-50].  It  is  interesting  to  note  that 
lithium  in  the  two  octahedral  positions  shows  distinct  6Li  chemical 
shifts.  The  difference  in  chemical  shift  between  these  two  octa¬ 
hedral  positions  can  be  explained  by  the  fact  that,  although  both 
16c  and  16d  lithium  sites  are  octahedrally  coordinated  by  oxygen, 
they  have  different  next-nearest  neighbors.  The  Li(16c)06  octa¬ 
hedron  is  in  a  corner-sharing  configuration  with  TiOg  octahedra, 
whereas  Li(16d)06  octahedra  are  in  corner-sharing  configurations 
with  vacant  and  Li(  1 6c)06  octahedra.  Such  a  difference  in  the  nature 
of  the  lithium  next-nearest  neighbors  is  likely  to  be  responsible 
for  the  observed  difference  in  chemical  shifts  between  the  16c  and 
16d  sites.  Indeed,  the  significant  effect  of  next-nearest  neighbors 
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has  been  documented  previously,  for  example,  for  23  Na  and  93  Nb 
chemical  shifts  [52,53]. 

Additionally,  of  all  the  NMR  spectra  collected  in  this  study,  that 
obtained  at  295  K  offers  the  best  resolution  of  the  three  resonance 
peaks.  The  integral  intensities  determined  from  this  spectrum  are 
62  ±3%,  24  ±3%,  and  14  ±2%  for  the  resonances  positioned  at 
0.65,  0.17,  and  -0.45  ppm,  respectively.  This  observation  provides 
further  support  for  the  assignment  presented  above  as  24%  corre¬ 
sponds  within  uncertainty  to  the  expected  proportion  of  Li  ions  in 
16d  sites,  whereas  inversing  the  assignment  of  the  two  octahedral 
sites  would  result  in  12%  of  the  Li  ions  in  16d  sites,  which,  in  turn, 
would  imply  an  unrealistically  large  number  of  vacancies  in  16d 
sites.  Finally,  as  explained  in  detail  below,  the  integral  intensity  of 
the  resonance  assigned  to  the  16c  site  increases  significantly  with 
temperature.  This  observation  is  consistent  with  our  assignment 
given  that  inversing  the  assignment  of  the  two  octahedral  sites 
would  signify  that  vacancies  at  16d  sites  are  progressively  filled 
by  lithium  ions  as  the  temperature  increases.  However,  this  phe¬ 
nomenon  would  go  directly  against  the  observations  of  Leonidov  et 
al.  [54],  who  reported  the  formation  of  vacancies  at  16d  sites  with 
increasing  temperature. 

Having  assigned  all  three  peaks,  the  variable-temperature  6 Li 
MAS  NMR  spectra  can  offer  further  information  on  the  possible 
lithium  diffusion  behavior  in  Li4Ti5  Ou  by  studying  the  temperature 
dependence  of  the  peak  positions.  Fig.  4a  shows  the  temperature 
dependence  of  the  6Li  chemical  shift  of  all  three  lithium  resonances. 
Two  general  phenomena  should  be  noted  to  interpret  the  NMR 
data.  Firstly,  thermal  expansion  of  a  material  can  lead  to  mono¬ 
tonic  changes  in  chemical  shift  due  to  changes  in  bond  lengths  and 
bond  angles  [55].  Secondly,  an  increasing  extent  of  lithium  hopping 
between  two  different  sites  leads  to  the  coalescence  of  the  corre¬ 
sponding  6Li  resonance  peaks.  The  main  characteristic  feature  of 
coalescence  is  the  change  in  isotropic  chemical  shift  observed  when 
two  peaks,  representing  two  different  chemical  sites,  move  toward 
each  other.  After  such  coalescence,  the  coalesced  peak  width  will 
narrow  with  increasing  temperature  due  to  motional  narrowing 
effects  [56]. 

As  can  be  seen  from  Fig.  4a,  all  three  peaks  show  shifts  towards 
lower  shielding  (or  decreased  ppm)  with  increasing  temperature. 
Therefore,  this  general  trend  can  be  attributed  to  the  thermal 
expansion  of  the  material  with  temperature.  Below  525  K,  the  8 a 
peak  shows  significant  change  in  its  isotropic  chemical  shift  and 
moves  towards  the  16d  peak  and  above  525  K  the  8 a  and  16d 
peaks  become  inseparable.  However,  this  may  not  represent  the 
coalescence  of  the  8a  and  16d  resonances  due  to  the  fact  that  the 
coalescence  should  be  followed  by  narrowing  of  the  coalesced  peak, 
but  all  the  peaks  show  an  increase  in  line  width  over  the  whole 
temperature  range.  These  observations  lead  us  to  conclude  that  the 
broad  peak  above  525 1<  is  not  the  result  of  the  coalescence  of  the 
8 a  and  1 6 d  peaks,  but  instead  is  due  to  the  broadening  of  the  8a  and 
16d  peaks  to  the  extent  that  they  strongly  overlap  and  thus  cannot 
be  clearly  resolved.  In  other  words,  there  is  no  significant  lithium 
exchange  between  the  8 a  and  16d  sites  which  could  lead  to  the 
coalescence  of  these  signals.  This  conclusion  is  supported  by  our 
molecular  simulations,  which  will  be  presented  later,  as  well  as  by 
previous  studies  reported  by  other  researchers  [50,54,57],  which 
indicate  that  lithium  ions  in  16d  sites  are  not  participating  in  any 
exchange  processes  at  these  temperatures. 

Deconvolution  of  the  8 a  and  1 6d  resonances  above  525  K,  would 
result  in  unacceptable  uncertainties;  therefore,  we  decided  to  fit 
the  broad  peak  with  a  single  resonance.  Above  525  K,  the  com¬ 
bined  8a/l  6d  peak  and  the  1 6c  peak  move  towards  higher  shielding 
(i.e.,  towards  negative  ppm).  At  the  same  time,  both  the  8a/16d 
and  16c  peaks  show  an  increase  in  line  width  with  temperature, 
an  indication  of  possible  coalescence.  However,  the  strong  overlap 
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Fig.  4.  Temperature  dependence  of  deconvoluted  parameters  (a)  isotropic  6Li  chem¬ 
ical  shift  and  (b)  integral  intensity  from  variable  temperature  6Li  MAS  NMR  spectra 
measured  at  21.1  T  with  spinning  frequency  of  101<Hz  for  Li4Ti50i2  (lines  are  to 
guide  the  eye).  The  peaks  8a  and  16d  are  inseparable  above  525  K,  hence  the  star 
symbol  represents  the  combination  of  both  8a  and  16d  peaks. 


between  the  8 a  and  16d  resonances  prevents  us  from  accurately 
determining  the  position  of  the  8 a  maximum.  As  a  result,  we  can¬ 
not  resolve  any  potential  coalescence  of  the  8  a  and  1 6c  resonances; 
however,  the  data  obtained  above  525 1<  do  not  rule  out  coalescence 
either.  Lithium  hopping  between  8a  and  16c  sites  is  expected  to 
take  place  based  on  our  molecular  dynamics  simulations  as  well 
as  previous  studies  [50,54,57].  Therefore,  in  order  to  explore  the 
lithium  exchange  between  the  8a  and  16c  sites  further,  we  ana¬ 
lyzed  the  temperature  dependence  of  the  integral  intensity  of  the 
6Li  resonances. 

Fig.  4b  shows  the  temperature  dependence  of  the  6Li  integral 
intensity  of  all  three  lithium  resonances.  The  integral  intensity  of 
the  16c  peak  increases  from  12%  at  300  K  to  about  40%  at  675  K 
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at  the  expense  of  the  8 a  peak,  which  decreases  drastically  with 
increasing  temperature.  This  result  clearly  indicates  a  progressive 
redistribution  of  the  lithium  ions  from  the  8a  site  to  the  16c  site 
as  the  temperature  increases.  Such  a  lithium  migration  should  ulti¬ 
mately  lead  to  the  coalescence  of  the  respective  6Li  peaks.  Indeed, 
the  complete  migration  of  lithium  from  8 a  to  16c  sites  would  lead 
to  a  defective  NaCl-type  structure  (space  group  Fm3m),  whereby 
lithium  ions  and  cation  vacancies  distribute  over  the  16c  sites. 

The  relocation  of  lithium  from  the  8 a  to  the  16c  site  with 
increasing  temperature  could  have  a  detrimental  effect  on  the  elec¬ 
trochemical  performance  of  Li4Ti50i2.  It  has  been  reported  that, 
as  a  result  of  electrochemical  lithium  insertion,  lithium  ions  move 
from  the  8a  to  16c  sites  [6,58]  and  consequently,  electrode  perfor¬ 
mance  is  dependent  on  the  number  of  available  16c  sites  in  the 
structure.  The  present  high-temperature  NMR  results  show  that  a 
similar  relocation  of  lithium  ions  from  8a  to  16c  sites  is  caused  by 
increasing  temperature.  This  relocation  could  adversely  affect  the 
number  of  free  1 6c  sites  available  for  electrochemical  lithium  inser¬ 
tion  and  thus  may  lead  to  poor  electrode  performance  of  Li4Ti50i2. 
However,  it  is  worthwhile  to  mention  that  this  effect  becomes 
significant  at  temperatures  well  above  the  normal  operating  tem¬ 
perature  of  lithium  batteries. 

Further,  it  is  interesting  to  note  that  there  is  no  significant 
change  in  the  integral  intensity  of  the  16d  peak  with  increasing 
temperature  (up  to  525 1<),  which  suggests  that  there  is  no  net 
lithium  migration  from  1 6 d  positions  to  other  sites  in  this  tempera¬ 
ture  range.  Such  a  lithium  migration  from  16d  to  16c  was  observed 
above  800  K  by  Leonidov  et  al.  [29]  leading  to  a  second  phase 
transformation  and  the  formation  of  an  ordered  NaCl-type  struc¬ 
ture.  Unfortunately,  due  to  technical  difficulties  in  building  a  very 
high  temperature  NMR  probe  (>700  K),  we  were  not  able  to  ver¬ 
ify  the  occurrence  of  this  second  phase  transformation  (where  all 
lithium  ions  in  1 6 d  sites  move  to  1 6c  sites).  Nonetheless,  the  present 
NMR  results  suggest  that  the  completion  of  a  phase  transforma¬ 
tion  of  Li4Ti50i2  from  a  spinel  to  a  defective  NaCl-type  structure 
should  occur  at  temperatures  above  700  K.  This  result  is  consistent 
with  the  previous  high-temperature  spectroscopic  studies  [28,29] 
of  Li4Ti50i2  discussed  in  Section  1. 

To  provide  further  evidence  for  the  lithium  diffusion  behavior 
determined  on  the  basis  of  6Li  MAS  NMR,  we  carried  out  two- 
dimensional  (2D)  6Li  exchange  spectroscopy  (EXSY)  [59].  Fig.  5 
shows  the  6Li  2D  EXSY  spectra  obtained  at  350  and  570  K  with  a 
mixing  time  (rmix)  of  1  s  and  a  spinning  speed  of  7  kHz  under  a  mag¬ 
netic  field  of  21.1  T.  In  2D  EXSY  spectra,  if  the  exchange  between 
any  two  sites  is  slow  relative  to  the  selected  mixing  time  then  peaks 
appear  only  on  the  diagonal  position.  On  the  other  hand,  if  the 
timescale  for  exchange  is  faster  than  rmix,  cross  peaks  appear  at 
the  co-ordinates  of  the  corresponding  diagonal  peaks  [59]. 

The  6Li  2D  EXSY  spectrum  obtained  at  350 1<  (see  Fig.  5a)  clearly 
shows  the  absence  of  any  cross  peaks.  The  detection  of  cross  peaks 
corresponding  to  exchange  between  the  8a  and  16c  resonances 
could  have  been  made  difficult  by  the  lithium  diffusion  pathways 
operating  in  this  material  at  this  temperature.  At  350  K,  more  than 
60%  of  the  lithium  ions  reside  in  8a  sites  whereas  the  16c  sites  only 
account  for  approximately  10%  of  the  lithium  ions.  Therefore,  most 
of  the  lithium  diffusion  involves  lithium  ions  diffusing  from  8a  to  8a 
sites.  However,  considering  the  crystal  structure  of  Li4Ti50i2  (see 
Fig.  1 ),  lithium  hopping  between  two  8 a  sites  has  to  take  places  via 
a  16c  site.  We  propose  that,  given  the  lithium  partition  between 
the  8 a  and  16c  sites,  lithium  occupancy  of  the  16c  site  is  expected 
to  be  transient  and  thus  would  not  be  detected  in  the  2D  EXSY 
experiments.  Indeed,  in  the  Li4Ti5012  structure,  a  16c  site  has  four 
nearest-neighbor  8 a  sites  at  1 .84  A.  In  addition,  based  on  the  lithium 
partition  between  the  8 a  and  16c  sites  derived  from  6Li  MAS  NMR, 
a  lithium  ion  in  a  16c  site  has,  on  average,  3.4  occupied  nearest- 
neighbor  8a  sites  at  350  K.  Such  a  close  proximity  between  lithium 


Fig.  5.  Two  dimensional  (2D)  6Li  exchange  spectroscopy  (EXSY)  spectra  of  Li4Ti50i2 
measured  with  a  21.1  T  magnetic  field  and  spinning  speed  of  7  kHz  at  350  K  (a)  and 
570  K  (b).  The  mixing  time  (rmiX)  is  kept  at  1  s  to  ensure  the  detecting  of  cross  peaks 
due  to  very  slow  jump  rates  (in  the  order  few  Hz).  The  single  pulse  6Li  MAS  NMR 
spectra  at  respective  temperature  are  plotted  on  the  top  of  the  2D  spectra  for  clear 
spectral  assignments.  The  arrows  indicate  the  expected  position  of  cross  peaks. 

ions  in  8a  and  16c  sites  makes  the  16c  site  energetically  costlier 
due  to  Li— Li  repulsive  electrostatic  interactions.  Such  a  diffusion 
behavior  has  already  been  proposed  by  Wagemaker  et  al.  [50],  who 
concluded,  based  on  the  results  of  ID  and  2D  7Li  MAS  NMR,  that 
lithium  ions  will  spend  negligible  time  in  16c  sites  while  diffusing 
between  8a  sites.  Additionally,  although  our  experiments  do  not 
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provide  further  direct  evidence,  it  is  worthwhile  to  mention  that 
our  results  are  consistent  with  the  nano-domain  model  proposed 
by  Wagemaker  et  al.  [50],  in  which  lithium  distributes  between 
separate  8a-  and  16c-rich  nano-domains. 

With  increasing  temperature,  the  proportion  of  lithium  ions  at 
8 a  sites  decreases  considerably  (see  Fig.  4b),  thereby  increasing  the 
probability  of  jumps  between  8a  and  16c  sites.  Any  such  jumps 
between  8a  and  16c  sites  can  be  detected  in  2D  EXSY  measure¬ 
ments.  Fig.  5b  shows  the  2D  EXSY  spectrum  measured  at  575  K. 
The  arrows  indicate  the  expected  positions  of  the  cross  peaks  that 
would  arise  from  any  diffusion  between  8 a  and  16c  sites.  Unfortu¬ 
nately,  identification  of  the  cross  peaks  is  difficult  due  to  the  severe 
overlap  between  the  6Li  resonances.  As  a  result,  the  low  spectral 
resolution  obtained  at  high  temperatures  prevented  us  from  clearly 
identifying  the  cross  peaks  and  calculating  the  lithium  jump  rate 
between  8a  and  1 6c  sites.  However,  it  should  be  noted  that  all  spec¬ 
tra  were  measured  at  magnetic  field  of  21 .1 T,  which  is  higher  than 
any  magnetic  field  used  in  studies  of  lithium  diffusion  in  this  mate¬ 
rial  published  to  date.  Therefore,  although  not  always  sufficient  for 
determining  all  the  properties  of  interest  with  high  accuracy,  we 
note  that  it  is  the  very  high  resolution  of  the  NMR  measurements 
presented  in  this  work  which  makes  this  study  possible. 

Based  on  our  conclusions  regarding  lithium  migration,  we  can 
attempt  to  explain  the  considerable  increase  in  the  conductiv¬ 
ity  of  Li4Ti50i2  at  high  temperatures  reported  in  the  literature 
[29].  The  integral  intensity  of  the  8 a  peak  decreases  from  62%  at 
room  temperature  to  48%  at  525  K  (see  Fig.  4b).  Such  migration 
leads  to  a  high  concentration  of  vacancies  at  8 a  sites  thereby  facil¬ 
itating  lithium  diffusion  through  16c^8a^l6c  pathways.  The 
16d  octahedral  site  shows  no  significant  lithium  diffusion  in  our 
NMR  spectra,  therefore  the  substantial  increase  in  the  conductivity 
reported  in  Li4Ti50i2  at  high  temperatures  [29]  is  caused  by  the 
increased  probability  of  lithium  diffusion  between  octahedral  16c 
sites  through  nearby  vacant  8a  sites.  In  order  to  cross-examine  the 
lithium  migration  and  diffusion  characteristics  observed  by  NMR  at 
high  temperatures,  we  carried  out  an  extensive  molecular  dynam¬ 
ics  study  of  lithium  diffusion  in  this  material. 

4.2.  Lithium  distribution  between  crystallographic  sites 

In  this  section,  we  present  results  obtained  from  molecular 
dynamics  simulations  of  Li4Ti50i2  carried  out  at  a  range  of  tem¬ 
peratures  in  order  to  complement  the  NMR  data  discussed  above. 
In  addition,  we  extended  our  study  to  consider  higher  lithium 
contents,  i.e.,  Li4+xTi50i2  with  0<x<3.  This  allows  for  a  direct 
comparison  with  published  experimental  data  of  Li  distribution 
between  crystallographic  sites  and  thus  constitutes  an  additional 
evaluation  of  the  potential  model.  Moreover,  as  the  lithium  content 
can  be  easily  varied  in  the  molecular  simulations,  it  allows  us  to 
investigate  how  the  conclusions  drawn  from  the  study  of  Li4Ti50i2 
extend  to  compounds  with  higher  lithium  contents.  Therefore,  four 
different  lithium  contents  were  considered,  namely,  Li4+xTi50i2 
with  x  equals  to  0,  1,  2,  or  3.  For  each  value  of  x,  a  series  of  MD 
simulations  was  carried  out  whereby  the  temperature  was  varied 
from  300  to  700  K  with  100K  intervals  to  reproduce  the  tempera¬ 
ture  range  used  for  the  NMR  measurements.  In  each  simulation,  The 
Ti3+/Ti4+  distribution  in  the  supercell  was  determined  randomly 
and  was  kept  fixed  during  the  MD  simulation.  For  each  simulation, 
the  lithium  distribution  between  the  three  crystallographic  posi¬ 
tions  (Fig.  1)  was  determined  as  follows:  (a)  the  average  titanium 
positions  were  computed  using  the  entire  run;  (b)  the  lithium  crys¬ 
tallographic  positions  were  calculated  from  the  average  titanium 
positions  and  the  known  crystal  symmetry;  (c)  radial  distribution 
functions  (RDFs)  were  constructed  using  the  instantaneous  lithium 
positions  relative  to  each  of  the  three  crystallographic  sites;  and 
(d)  an  integration  of  each  RDF  plot  to  its  first  minimum  yielded  the 


Table  2 

Population  of  each  crystallographic  site  as  a  function  of  temperature  and  lithium 
content  as  determined  from  the  molecular  dynamics  simulations. 


Li  content,  x 

300  K 

400  K 

500  K 

600  K 

700  K 

8a 

0 

74 

74 

74 

73 

73 

1 

50 

49 

49 

49 

49 

2 

30 

29 

28 

28 

28 

3 

13 

12 

12 

11 

12 

16c 

0 

1 

1 

3 

4 

4 

1 

31 

31 

31 

32 

32 

2 

54 

54 

56 

56 

57 

3 

73 

74 

74 

75 

75 

16d 

0 

25 

25 

24 

24 

23 

1 

20 

20 

20 

20 

20 

2 

17 

17 

17 

17 

17 

3 

14 

14 

14 

14 

14 

occupancy  of  the  three  sites.  Table  2  shows  the  lithium  distribu¬ 
tions  at  all  temperatures  and  all  Li  contents  and  Fig.  6  compares 
the  lithium  distribution  at  300  K  with  the  experimental  data  of 
Scharner  et  al.  [24],  Wagemaker  et  al.  [25],  and  Kataoka  et  al.  [42]. 
The  MD  simulations  predict  a  full  occupancy  of  the  16d  sites  based 
on  the  stoichiometry  at  all  temperatures  in  agreement  with  the 
experimental  data  of  Wagemaker  et  al.  and  Kataoka  et  al.  as  well 
as  our  NMR  data.  The  potential  model  also  shows  a  Li  distribu¬ 
tion  between  the  8 a  and  16c  sites  in  excellent  agreement,  for  all 
lithium  contents,  with  the  room  temperature  experimental  data  of 
Scharner  et  al.  and  Kataoka  et  al.,  with  only  a  slight  overestima¬ 
tion  of  the  occupancy  of  the  8a  site  at  higher  Li  contents.  The  data 
of  Wagemaker  et  al.  is  shifted  to  higher  Li  contents  but  follow  the 
same  trend  as  the  XRD  and  MD  results. 

Turning  now  to  the  temperature  dependence,  our  MD  simu¬ 
lations  show  an  increase  in  the  occupancy  of  the  16c  site  to  the 
detriment  of  the  8 a  site  with  increasing  temperature  (Table  2). 
Although  this  trend  is  in  agreement  with  the  NMR  data  obtained  in 
this  work,  it  is  not  as  pronounced  as  observed  experimentally.  As 
mentioned  previously,  Raman  [29]  and  infrared  [28]  spectroscopy 
experiments  also  indicate  the  migration  of  lithium  ions  from  the  8 a 
to  the  16c  site  with  increasing  temperature.  To  conclude,  it  is  not 
clear  why  the  lithium  content  dependence  of  the  Li  distribution  is 
well  reproduced  whereas  the  temperature  dependence  is  far  too 
weak  compared  to  experiment.  This  result  might  suggest  that  the 


Fig.  6.  Lithium  distribution  between  the  8 a  and  16c  sites  as  a  function  of  Li  content 
as  obtained  from  molecular  dynamics  simulations  (this  work),  neutron  diffraction 
(ND  [25]),  and  X-ray  diffraction  data  (XRD1  -  Kataoka  etal.  [42]  and  XRD2  -  Scharner 
et  al.  [24]). 
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strength  of  the  Li— Li  repulsion  relative  to  the  Li-0  potential  is  well 
reproduced  by  the  potential  model  but  that  the  Li-0  interaction  is 
too  attractive  with  respect  to  the  thermal  energy. 


4.3.  Lithium  diffusion  between  crystallographic  sites 


The  same  series  of  MD  simulations  was  used  to  extract  informa¬ 
tion  about  the  rate  and  mechanism  of  lithium  diffusion  for  each  of 
the  four  lithium  contents.  As  will  be  demonstrated  in  this  subsec¬ 
tion,  the  rates  of  exchange  between  the  different  crystallographic 
sites  span  a  large  range  of  timescales.  Therefore,  two  approaches 
were  used  to  determine  the  rates  of  lithium  exchange  as  a  function 
of  temperature  and  thus  extract  activation  energies  to  be  compared 
with  experimental  data. 

The  first  approach  consists  of  calculating  the  population  cor¬ 
relation  function,  C(t),  for  lithium  hopping  between  pairs  of 
crystallographic  sites.  C(t)  is  defined  as  [60]: 


(hA(0)hB(t)) 

(hA( 0)) 


(3) 


where  hfx)  is  1  if  the  lithium  ion  is  in  site  i  at  time  x  and  0  oth¬ 
erwise.  Correlation  was  followed  for  20  ps  and  C(t)  was  obtained 
by  averaging  the  contribution  of  all  lithium  ions.  As  explained  by 
Chandler  [60],  this  correlation  function  approaches  its  asymptotic 
value  exponentially: 


1  -  exp 


(4) 


where  r is  the  inverse  of  the  sum  of  the  forward  (kAB)  and 
backward  (kBA)  reaction  rates.  In  addition,  if  r rxn  is  much  larger 
than  the  molecular  time  scale  ( r moi ),  C(  t)  is  linear  in  the  time  regime 
rmoi  <t<Trxn  and  its  gradient  is  equal  to  kAB. 

The  second  approach  consists  of  computing  the  residence  time 
correlation  function,  R(t): 


where  Nt  is  the  number  of  lithium  ions  in  the  crystallographic  site  of 
interest  at  time  t  and  Oft)  is  the  Heaviside  function,  which  is  1  if  the 
ith  lithium  ion  is  in  the  crystallographic  site  of  interest  at  time  t  and 
0  otherwise.  The  crystallographic  sites  were  considered  to  extend  to 
the  first  minima  of  the  radial  distribution  functions  described  in  the 
previous  subsection.  As  it  is  usually  done  with  such  an  analysis  [61  ], 
a  lithium  ion  was  counted  as  having  left  its  crystallographic  site  if  it 
has  done  so  for  any  continuous  period  of  time  of  at  least  t*.  This  is  to 
allow  for  lithium  ions  that  temporarily  leave  the  crystallographic 
site  and  come  back  without  visiting  another  site  for  a  significant 
amount  of  time  to  be  counted  in  the  residence  time.  A  value  of 
0.5  ps  for  t*  was  employed  in  this  work.  The  residence  time  and 
thus  the  exchange  rate,  k,  can  then  be  obtained  by  integration  of 
<R(t)): 


0.05 


0  5  10  15  20 

Time  (ps) 


1/T  (K'1) 


poo 

Ti=r  =  J  (m)dt  (6) 

The  first  approach  is  appropriate  when  r moi «:  r rxn  whereas  the 
second  approach  converges  faster  the  faster  the  exchange  rate. 

The  population  correlation  functions  were  calculated  for  lithium 
exchange  between  8a  and  16d,  16d  and  16c,  and  16c  and  8a.  These 
functions  revealed  that  there  are  no  direct  exchanges  between  the 
8 a  and  16d  sites.  Indeed,  inspection  of  the  crystal  structure  shows 
that  the  1 6c  site  is  closer  to  1 6d  than  to  8a  (see  Fig.  1 )  and  the  popu¬ 
lation  correlation  function  therefore  confirms  the  expectation  that 
the  lithium  ions  hop  first  to  1 6c  from  1 6 d.  As  expected,  the  1 6d-l  6c 
correlation  function  shows  a  linear  behavior  at  all  temperatures 


Fig.  7.  Population  correlation  function  for  the  16d-16c  Li  exchange  as  a  function  of 
temperature  (top);  rate  of  Li  exchange  between  the  16d  and  16c  sites  as  a  function 
of  temperature  for  each  Li  content  (middle);  activation  energy  for  the  16d-16c  Li 
exchange  as  a  function  of  Li  content  (bottom). 

(see  Fig.  7  for  an  example),  whereby  the  gradient  of  each  straight 
line  yields  the  forward  rate  constant.  A  plot  of  the  logarithm  of 
the  rate  constant  as  a  function  of  the  inverse  temperature  (Fig.  7) 
gives  the  activation  energy  for  each  Li  content.  Interestingly,  Fig.  7 
also  shows  that  the  activation  energy  for  lithium  hopping  between 
the  16d  and  16c  sites  is  a  linear  function  of  the  lithium  content 
with  activation  energies  from  approximately  43kJmol-1  for  x  =  0 
to  63  kj  mol-1  for  x  =  3. 
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Fig.  8.  Li  residence  time  in  the  1 6c  site  for  x  =  3  at  all  temperatures  (top)  and  rate  of 
exchange  as  a  function  of  temperature  for  all  Li  contents  (bottom). 

The  exchange  between  the  1 6c  and  8 a  sites  was  fast  with  no  dis¬ 
tinct  linear  regime.  Therefore,  the  lithium  residence  time  in  each 
of  the  8a  and  16c  sites  was  calculated  using  the  second  approach. 
The  residence  time  correlation  function  of  lithium  in  the  1 6c  site  for 
x  =  3  at  all  temperatures  is  shown  in  Fig.  8  as  an  example.  The  rate 
of  exchange  is  then  obtained  from  Eq.  (6).  The  activation  energies 
for  hopping  from  both  sites  initially  increase,  reach  a  maximum 
and  then  decrease  again  with  increasing  Li  content.  However,  the 
maximum  activation  energies  are  obtained  at  different  Li  contents. 
The  calculated  activation  barriers  are  of  the  order  of  5-10  kj  mol-1 
and  are  therefore  much  lower  than  the  activation  energy  values 
of  approximately  30  kjmol-1  derived  by  Wagemaker  et  al.  [27] 
from  NMR  correlation  times  of  lithium  hopping  on  the  8a  and  16c 
sublattices.  The  activation  barriers  obtained  by  Wagemaker  et  al. 
were  for  mixed  8a/16c  occupation  (i.e.,  x  =  0.3,  1,  and  2)  and  no 
lithium  diffusion  was  observed  for  the  structural  end  members, 
x  =  0  and  3.  This  observation  is  consistent  with  the  ultraslow  Li  dif¬ 
fusion  reported  by  Wilkening  et  al.  [30,31  ]  for  x  =  0.  In  our  molecular 
dynamics  calculations,  it  is  likely  that  a  hop  out  of  either  site  is 
often  followed  by  a  backward  hop  to  the  same  site  thus  leading  to 
no  net  lithium  diffusion.  To  address  this  issue,  we  calculated  the 
overall  lithium  diffusion  coefficient  as  a  function  of  temperature  to 
extract  mean  activation  energies  for  lithium  diffusion  whereby  the 
lithium  diffusion  coefficients  were  obtained  from  the  gradients  of 
the  mean  square  displacement  functions.  Table  3  shows  the  overall 


Table  3 

Activation  energies  (kjmol-1)  for  lithium  hopping  from  the  8a  and  16c  sites  in 
Li4+xTi50i2  as  determined  from  residence  time  correlation  functions  and  activation 
energies  for  overall  lithium  diffusion  as  obtained  from  mean  square  displacement 
functions. 


Site 

LUTisO^ 

Li5Ti50i2 

LieTisO^ 

LiyTisO^ 

8a 

4.2 

7.0 

9.3 

8.5 

16c 

5.2 

8.6 

7.3 

5.3 

Overall 

26.8 

21.5 

20.3 

19.6 

activation  energies  as  a  function  of  Li  content.  The  activation  ener¬ 
gies  obtained  from  the  MSD  functions  are  much  higher  than  those 
computed  from  the  residence  time  correlation  functions,  which 
confirm  the  hypothesis  stated  above.  The  activation  energies  for 
the  mixed  occupation  values  are  lower  than  those  obtained  exper¬ 
imentally  by  Wagemaker  et  al.  [27]  (~30  kj  mol-1 )  and  Wilkening 
et  al.  [30]  (~40  kjmol-1).  In  addition,  the  activation  energies  for 
the  structural  end  members  are  much  lower  than  those  derived 
experimentally  (70-90 kjmol-1)  [31].  However,  the  weak  depen¬ 
dence  of  the  calculated  activation  energies  on  the  Li  content  for 
the  mixed  occupation  compositions  agrees  with  the  experimental 
observation  of  Wagemaker  et  al.  [27]. 

In  conclusion,  the  molecular  dynamics  simulations  show  that 
the  free  energy  barrier  for  lithium  hopping  between  the  8a  and 
16c  sites  is  much  lower  than  that  between  the  16d  and  16c  sites 
and  that  therefore  the  onset  of  the  first  transfer  should  occur  at 
a  lower  temperature.  This  is  in  agreement  with  the  NMR  findings 
and,  although  the  motion  of  the  Li  ions  in  16d  positions  was  not 
observed  in  the  NMR  experiments,  the  work  of  Leonidov  et  al.  [29] 
indicates  that  the  1 6 d  Li  ions  begin  to  move  at  temperatures  higher 
that  those  probed  in  this  work. 

5.  Conclusions 

Highly  crystalline  Li4Ti50i2  has  been  synthesized  by  a  modified 
non-alkoxide-based  sol-gel  method,  where  the  final  calcination 
is  performed  at  a  lower  temperature  and  for  a  shorter  period  of 
time  than  in  previously  reported  sol-gel  methods.  This  sol-gel  syn¬ 
thesis  method  constitutes  a  simple  and  cost  efficient  preparation 
method  suitable  for  bulk  industrial  preparations.  High  temperature 
6Li  MAS  NMR  measurements  of  the  synthesized  Li4Ti50i2  sample 
reveal  the  migration  of  lithium  ions  from  the  tetrahedral  8a  sites 
to  the  octahedral  16c  sites  leading  to  a  defective  NaCl-type  struc¬ 
ture.  This  lithium  relocation  could  restrict  the  number  of  free  16c 
sites  available  for  external  lithium  insertion  and  may  lead  to  poor 
electrode  performance  of  Li4Ti50i2.  Due  to  this  phase  transfor¬ 
mation,  the  operating  temperature  of  any  Li4Ti50i2  based  lithium 
battery  should  be  kept  below  500  K.  Interestingly,  the  lithium  ions 
in  octahedral  16d  sites  show  no  significant  diffusion  process  in  the 
experimental  temperature  range.  The  change  in  lithium  population 
observed  between  the  8 a  and  1 6c  peaks  also  suggests  lithium  diffu¬ 
sion  through  16c^  8a  ->  16c  pathways,  which  is  facilitated  by  the 
large  number  of  vacancies  formed  at  8a  sites.  Therefore,  the  large 
increase  in  conductivity  observed  for  Li4Ti50i2  at  high  tempera¬ 
tures  (>600 1<)  is  mainly  due  to  the  high  lithium  diffusion  between 
16c  sites  through  vacant  8a  sites. 

Our  potential  shell  model  gives  a  good  structural  description 
of  the  spinel  oxides  Li4+xTi50i2  (0<x<3)  and  an  excellent  agree¬ 
ment  with  the  experimentally-derived  Li  distributions  between  the 
8a  and  16c  sites  as  a  function  of  Li  content.  However,  the  depth  of 
the  free  energy  minima  of  either  or  both  of  the  8a  and  16c  posi¬ 
tions  is  probably  too  shallow,  which  leads  to  a  higher  Li  diffusion 
and  is  also  likely  to  be  the  cause  of  the  smaller  temperature  depen¬ 
dence  of  the  Li  distribution  between  the  8a  and  16c  sites.  Despite 
the  discrepancies  in  the  absolute  values  of  the  activation  energy  for 
lithium  diffusion,  important  information  emerged  from  the  molec- 
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ular  dynamics  simulations  concerning  the  relative  dynamics  of 
exchange  between  the  different  lithium  sites  in  support  of  our  high 
temperature  6Li  MAS  NMR  data.  At  first,  the  exchange  between  the 
8 a  and  16c  sites  is  relatively  faster  whereas  that  between  the  16d 
and  16c  sites  is  slow  and  that  there  is  no  direct  exchange  between 
the  16d  and  8 a  sites.  Additionally,  the  occupancy  of  the  16d  sites 
remains  at  nearly  the  stoichiometric  level  at  all  temperatures  stud¬ 
ied  as  it  is  found  in  the  NMR  analysis  of  the  relative  peak  intensities. 
Finally,  this  study  demonstrates  the  importance  and  effectiveness 
of  combining  information  obtained  from  NMR  measurements  and 
MD  simulations  to  draw  a  comprehensive  picture  of  the  micro¬ 
scopic  lithium  diffusion  processes  in  electrode  materials. 
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